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Quantum Random Access Codes using Single d-level Systems 
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Random access codes (RACs) are used by a party to despite limited communication access an 
arbitrary subset of information held by another party. Quantum resources are known to enable 
RACs that break classical limitations. Here, we study quantum and classical RACs with high-level 
communication. We derive average performances of classical RACs and present families of high-level 
quantum RACs. Our results show that high-level quantum systems can significantly increase the 
advantage of quantum RACs over the classical counterparts. We demonstrate our findings in an 
experimental realization of a quantum RAC with four-level communication. 


Introduction. — Quantum information arises from the 
realization that information can be encoded in physi¬ 
cal systems subject to the laws of quantum mechanics, 
and thus allowing for the use of non-classical phenom¬ 
ena such as superpositions and entanglement as quan¬ 
tum resources for information processing. Today, it is 
well-known that quantum resources can enhance informa¬ 
tion processing beyond classical limitations in cryptogra¬ 
phy [1-4], communication complexity [5-7] and numerous 
other information or computation problems [8-11]. 

A broad category of communication tasks that are use¬ 
ful for a wide variety of applications is random access 
codes (RACs). In a RAC, one party Alice holds a set of 
information and another party Bob aims to access an ar¬ 
bitrary subset of the information held by Alice. However, 
the communication between the parties is restricted. We 
can denote a RAC as n 1, meaning that Alice has n 
bits of information while Bob is interested in one of Al¬ 
ice’s bits. The task is for Alice to encode her n bits into 
only one bit, which she sends to Bob, such that Bob can 
find the value of any of her bits with a high probability 
P- 

In the analog quantum random access code (QRAC), 
we allow Alice to encode her n bits into a single two-level 
quantum system that she sends to Bob. Depending on 
which of Alice’s bits Bob wishes to find, he performs a 
suitable measurement on the system and with high prob¬ 
ability recovers the true value of the particular bit. One 
should distinguish between the following two scenarios: 
Bob recovers the bit of is interest with (i) a worst case 
probability p, and (ii) an average probability p. 

QRACs were initially introduced in studies of quan¬ 
tum finite automata [12-14]. Nevertheless, QRACs were 
shown useful for various other applications including lo¬ 
cally decodable codes [15], network coding [16], reduction 
of communication complexity [17], semi-device indepen¬ 
dent random number expansion [18] and semi-device in¬ 
dependent key distribution [19]. Additionally, QRACs 
are interesting for foundational studies of quantum and 
no-signaling resources [20, 21]. Furthermore, random ac¬ 
cess coding can be realized by means of entanglement 
and classical communication [22] and such entanglement- 
assisted RACs have been applied to enhance strategies in 


games [23]. 

In [12], an example of a 2 1 QRAC was presented 

which provides an advantage over the corresponding clas¬ 
sical RAC. Classically, one can achieve an average success 
probability = 3/4 which can be obtained as follows. 
Let Alice send Bob the value of her first bit every time. 
Whenever Bob is interested in this bit he recovers it with 
certainty and whenever he is interested in Alice’s second 
bit he is forced to guess, succeeding in half the cases, 
leading to p^ = 3/4. However, in the QRAC, Alice en¬ 
codes her two bits xq, xi € {0,1} into a two-level system 

IV’xo.i) = ^(|0) + ^ ((-!)"“+*(-l)"0|l))- If Bob 
is interested in xg he performs a measurement in the ba¬ 
sis (1)1)) (1)~1)} and if he is interested in xi 

he measures in the basis (l,i), ^ (1,—i)}. Since 
the encoding states are symmetrically distributed on the 
equator of the Bloch sphere such that they form the ver¬ 
tices of a square, the probabilities of all outcomes are 
the same, namely p^ = \ (l-|-l/-\/2) « 0.854. Thus, 
the QRAC outperforms the RAC, and we quantify the 
advantage by the ratio p^ jp'^ « 1.138. 

As mentioned in [12], the above example was extended 
by Chuang to a 3 1 QRAC using eight encoding 

states symmetrically distributed over the surface of the 
Bloch sphere forming the vertices of a cube, and adding 
the third measurement basis {(1,0), (0,1)}. Due to the 
symmetry, the success probability is always = 1/2(1-|- 
l/v/l) « 0.789 whereas classically, one can still achieve 
p^ = 3/4 with a somewhat more sophisticated strategy 
than in the 2 —>■ 1 case. Hence, the advantage of this 
QRAC over the RAC is p'^/p^ « 1.052. 

However, RACs (QRACs) can also be considered when 
Alice and Bob do not use bits (qubits) but a classical 
(quantum) d-level system. To the authors’ knowledge, 
little attention has been directed to such generalizations 
but initial steps have been taken in [17, 24]. In this letter, 
we study RACs and QRACs with high-level communica¬ 
tion, and we infer the notation 1 as an abbre¬ 

viation of the scenario where Alice has n d-levels i.e. a 
string X = xg...Xn-i where Xi G (0,..., d — 1}, which she 
encodes into a classical (quantum) d-level, which is sent 
to Bob who should be able to recover the Xj of his interest 
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with an average (or worst case) probability of Pn,d- 

Classical RAGs .— Let us begin by considering clas¬ 
sical RACs on the form ^ 1. Any deterministic 
classical strategy for the RAC is composed of an encod¬ 
ing function used by Alice to encode her string into a 
d-level, and a decoding function used by Bob to decode 
the obtained d-level with respect to Xj. Proving the op¬ 
timality of a classical strategy is difficult for the general 
nC) 1 setting. However, intuition strongly suggests 
that an optimal strategy is for Alice to use majority en¬ 
coding i.e. Alice counts the number of times each of the 
values {0, 1} appears in her string x and outputs 

the d-level that occurs most frequently. In case of a tie, 
she can output either d-level. Bob does identity decoding 
and thus outputs whatever he receives from Alice. This 
strategy was proven optimal for d = 2 in [25]. 

Now, we find the average classical success probability 
Pnd arbitrary n and d. Let all positive integer parti¬ 
tions of the number n be stored into the set X and let Xj 
denote the j’th partition. The total number of partitions 
is |A|. In total, there are d” strings that Alice could 
be holding, and for each string x we make a frequency 
list i.e. a table where we count the number of times any 
same d-level occurs in x. Thus, for any given x there 
exists a unique j such that x associates to the partition 
Xj-. a: ~ Xj. Consider the number of strings Nx- such 
that X ^ Xj. We compute Nxj combinatorially by 

n(<i-m+i)(i) 

\lCX, / rn=l 

where fci,..., | are the elements of Xj, Xj is the set of 

distinct elements occurring in Xj and Cxj (1) is the num¬ 
ber of times I appears in Xj. Given an ’’empty” string 
X with n positions we need to assign one of d possible 
values to the first ki positions in x. Similarly, we can 
assign any of the remaining d — I values to the next ^2 
positions in x. This process is continued until we have 
assigned one of the remaining d — \Xj \ -|- 1 values to the 
final k\x-\ positions in x. The number of such strings is 
given by the third factor in (1). However, if we arbitrar¬ 
ily permute x into x' , it still holds that x' ^ Xj, and 
therefore we multiply by the first factor in (I) which is 
the number of permutations of x. Nevertheless, if it hap¬ 
pens that ki occurs more than once in Xj, some strings 
will be counted Cxj(ki)\ times. By multiplication with 
the second factor in (1), each string is counted only once 
and thus (1) is the number of strings x ^ Xj. 

Using majority encoding, Alice will compute the max¬ 
imum frequency ma,x{Xj} and output the associated di- 
level to Bob who’s identity decoding has a probability of 
max{Aj}/n to return the correct value of the d-level of 


his interest. Therefore, the average success probability is 
I 

max{X,} (2) 

i=i 

Evaluating this expression for the general (n, d)-setting 
is cumbersome. Therefore we have computed two special 
cases by fixing n = 2, 3 

Observe that in the limit of d —)■ oo the classical success 
probability converges to I/n which corresponds to the 
trivial strategy where Bob simply guesses which of the n 
d-levels he was given. 

Quantum advantage of QRACs in [24 -]. — Let us now 
shortly apply our formula (2) to investigate the advan¬ 
tage of some known high-level QRACs. In [24], a set of 
QRACs on the form (d -I- —>■ I for d = 2,3,4, 5, 7, 8 

were constructed from states maximizing the discrete 
Wigner function. The authors compute the average suc¬ 
cess probability of their QRACs but are unable to com¬ 
pare these with the success probability of the correspond¬ 
ing RAC. Having found the success probability of the 
RACs, we can implement our results with the QRACs 
of [24] to quantify the quantum advantage, see table 1. 
In the case of d = 2 the result reproduces the known 


d 

in [24] 


P^/P^ 

2 

0.789 

0.75 

1.052 

3 

0.637 

0.593 

1.075 

4 

0.5424 

0.4961 

1.0933 

5 

0.4700 

0.4291 

1.0953 

7 

0.3720 

0.3420 

1.0876 

8 

0.3372 

0.3118 

1.0815 


TABLE I. Comparison of quantum and classical success prob¬ 
abilities for QRACs in [24]. 

3(2) 

—7^ 1 QRAC of Chuang. However, we see that the 
QRACs with n = d -I- 1 > 2 can provide an advantage 
larger than the previously known case of d = 2. 

High-level QRACs with n = 2 .— We will now con¬ 
struct a family of QRACs on the form —)■ I i.e. Alice 

has Xo,Xi G {0,..., d — 1} which she encodes into a quan¬ 
tum d-level system sent to Bob. For this purpose we 
will use two mutually unbiased bases (MUBs) [26, 27], 
namely the computational basis {[?)}; s-nd the fourier 
basis Je;) = where oj = e^ . Let us 

assume that Alice encodes the string xgxi = 00 into the 
state l^/’oo) = where IV 2 ,d = ^^2 -h is the d- 

dependent normalization. We now use the unitary oper¬ 
ators X = X)fc=o 1^ + l)(^l Z = to 
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define the general encoding state \ 4 >xoxi) as 

( 4 ) 

Bob will perform a measurement in the basis {\l)}i if he is 
interested in xq and in the basis {le;)}/ if he is interested 
in xi- We proceed by deriving the success probability 
p^d of this family of QRACs. 

Firstly, we expand the encoding state (4) onto the form 

(ko) + ^ Ei=d + 2^o))- Then, if 

Bob performs a computational basis measurement and 
his outcome is denoted by /, the probability distribution 
Pq of Bob’s outcomes is 



Po{i) = \mxox,)f 



Sl,Xo + 


^Xi(l-Xo) 

Vd 


2 


( 5 ) 


On the other hand, assume that Bob is interested in xi, 
then he will measure in the fourier basis giving rise to a 
probability distribution Pi, 


Pi(0 = |(e,|V’.ox,)P 




, — lxo 


Vd 


+ UJ~ 


Oxi,l 


( 6 ) 


Bob’s probability of success amounts to his measurement 
outcome indeed being the correct value of the Xj he is 
interested in. This corresponds to the outcome I = xq in 
(5) and the outcome ^ = xi in (6). Evaluating both these 
probabilities expressions will yield 

For any d the QRAC is non-trivial since > 1/2. In 

the case of d = 2 the well known result of P 22 ~ 0.854 is 
reproduced. Furthermore, since the success probability of 
Bob only depends on d, the proposed family of QRACs 
perform equally well both on average and in the worst 
case setting i.e. the probability to recover xq is always 
the same as the probability to recover xi, independent of 
the encoding or the measurement. 

The quantum success probability (7) can be compared 
to the classical success probability (3) . Evidently, it holds 
that P 2 ^ > P 2 ^ for all d. We quantify the advantage of 
the QRAC over the RAC by the ratio and in 

figure (1), we represent the quantum advantage by the 
red dots. The optimal quantum advantage is found at 
d = 6 at which p^g/p^g = I ~ 1-207. This 

is constitutes a significant improvement over the d = 2 
case where the quantum advantage is only about 1.138. 
In terms of absolute numbers, the difference between the 
quantum and the classical success probability is maxi¬ 
mized for d = 4 where P 24 — p^U = 1/8. 

Experimental realization.-Now, we will now present an 
experimental realization of the 2*^"^) —>■ 1 QRAC. In our 
experiment, the physical systems are defined by single 


FIG. 1. The ratio between the quantum and classical success 
probabilities as a function of d for the 2^“*^ —^ 1 (red) and 
1 (blue) QRACs. 


photons in polarization and path setup. The information 
is encoded in four basis states: |1) = |id, a), |2) = |F, a), 
|3) = \H,b) and |0) = \V,b), where (H) and (F) are 
horizontal and vertical polarization photonic modes re¬ 
spectively, and (a and b) are two spatial photonic modes 
of single photons. Any ququart state can be written as 
a|id, a) + j3\V,a) P^\H, b) -|-(5|F, 6). 

The single photons are generated from a heralded 
single photon source using a 390 nm pulsed laser pump¬ 
ing a 2 mm BBO crystal and creating twin photons 
through a spontaneous parametric down-conversion 
(SPDC) process. The idler photon is used as trigger 
and detected by a single photon detector Dt- To 
exactly define the spatial and spectral properties of the 
signal photon, the emitted photon modes are coupled 
into a single mode fiber (SMF) and passed through 
a 3 nm narrowband interference filter (F). We have 
characterized our heralded single photon source and 
we have found that the ratio between the coincidence 
photon counting due to single- and multi-photon pair 
emissions is below 0.1%. To prepare the initial state 
of the photon in \H), the signal photon is passing 
through a polarizer oriented to horizontal polarization 
direction. Alice’s photonic states were prepared by 
three suitably oriented half-wave plates HWP(di), 
HWP( 02 ) and HWP( 03 ), polarization beam splitter 
(PBS), and a setting of a phase shift PS((/) (see figure 
2), cos{ 29 i)cos{ 292 )\H,a) + cos{ 29 i)sin{ 292 )\V,a) + 

\sin{29i)sin{29z)\H, b) — sin{ 29 i)cos{ 292 ,)\V,b)]. 
Thus, by adjusting the HWP orientation angles 9i, Alice 
could produce any of the 16 required states I'i/xoxi) with 
xo, xi G {0,1,2,3}, see Table II for all settings. 

Bob chooses between two measurement settings, either 
in the basis of eigenstates of Z (computational basis) 
or X (Fourier basis). For the computational basis, the 
measurement consists in each of two spatial modes (a) 
and (b), of PBS and two single photon detectors Dzi 
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FIG. 2. Experimental set-up for the realization of four-level 
encoding QRAC. Alice’s states are encoded in single photon 
horizontal and vertical polarizations, and in two spatial modes 
and prepared by suitable orientation of three half-wave plates 
(HWP) at angles di (with i = 1,2,3), PBS and phase shift 
PS When Bob measures in the computational basis, the 
measurement consists, in each of two spacial modes (a) and 
(b), of PBS and two single photon detectors Dzi (i = 1, ...4). 
Bob’s measurement in the Fourier basis consists of an inter¬ 
ference at BS and at each of the outputs of the BS, of HWP 
oriented at 22.5, PBS, and two single photon detectors Dxi 
(i = l,...4) . 

(i = 1, ...4). The choice to measure in a particular basis 
is implemented by moving the mirrors (M) in and out 
with help of pico motors translation stages. The Fourier 
basis measurement consists of an interference at BS and 
at each of the outputs of the BS, of HWP oriented at 22.5, 
PBS, and two single photon detectors (Dxi (i = 1, -.-d). 

Our single-photon detectors (D^, i = 1,2,3,4) for the 
two measurements were Silicon avalanche photodiodes 
with detection efficiency r]d = 0.55, the dead time 50 
ns, and dark counts rate Rd — 400 Hz. All coincidence 
counts between the signal and idler photons are registered 
using an Multi-channel coincidence logic with a time win¬ 
dow of 1.7 ns. The number of detected photons was ap¬ 
proximately 2500 per second and the total time used for 
each experimental settings was 10 s. 

The probabilities needed for the QRAC are obtained 
from the number of detections in the single photon de¬ 
tectors. All experimental results are presented in Ta¬ 
ble H. The results are in very good agreement with the 
predictions of quantum mechanics, namely for an ideal 
experiment = 0.75. Our results yield an average 
quantum success probability 0.754±0.038 which violates 
the classical bound P'^ = 0.625. The errors come from 
Poissonian counting statistics and systematic errors. The 
main sources of systematic errors are the slight intrinsic 
imperfections of the PBSs and HWPs. 

High-level QRACs n = 3.— Let us consider a sec¬ 
ond family of QRACs, on the form 3^“^^ —>■ 1 i.e. Al- 



Ol 

02 

03 

4> 

Pexp 

Piip 


12.05 

22.5 

-9.22 

0 

0.747 ±0.036 

0.752 ±0.038 


12.05 

22.5 

-9.22 

TV 

0.748 ± 0.036 

0.774 ±0.038 

l/'02 

12.05 

-22.5 

9.22 

0 

0.749 ± 0.036 

0.752 ±0.039 


12.05 

-22.5 

9.22 

TV 

0.753 ± 0.036 

0.713 ±0.039 

i>io 

12.05 

22.5 

-35.78 

0 

0.766 ± 0.036 

0.763 ± 0.038 

i/'ll 

12.05 

22.5 

-35.78 

TV 

0.767 ±0.036 

0.750 ± 0.038 

i>i2 

12.05 

-22.5 

35.78 

0 

0.764 ± 0.036 

0.751 ±0.039 

'IplS 

12.05 

-22.5 

35.78 

TV 

0.764 ± 0.036 

0.748 ± 0.039 


-22.5 

32.95 

9.22 

0 

0.755 ± 0.036 

0.764 ±0.038 

i’21 

-22.5 

32.95 

9.22 

TV 

0.761 ±0.036 

0.754 ±0.038 

‘tp22 

22.5 

32.95 

-9.22 

0 

0.786 ± 0.036 

0.775 ±0.039 

'!/’23 

22.5 

32.95 

-9.22 

TV 

0.786 ± 0.036 

0.728 ± 0.039 

i>30 

-22.5 

32.95 

35.78 

0 

0.742 ± 0.037 

0.772 ±0.038 

l/'Sl 

-22-5 

32.95 

35.78 

TV 

0.743 ± 0.037 

0.766 ± 0.038 

i>32 

22.5 

32.95 

-35.78 

0 

0.732 ± 0.037 

0.749 ± 0.039 

i>33 

22.5 

32.95 

-35.78 

TV 

0.733 ± 0.037 

0.703 ± 0.039 


TABLE II. The orientation of the three half wave plates 6i 
(with i = 1,2, 3) and the phase shift PS ((/>) for the 16 quan¬ 
tum states The quantum success probabilities P^p and 
Pexp for measurements in the computational and Fourier basis 
respectively. 


ice encodes xo,xi,X 2 S {0, 1} into a quantum 

d-level system sent to Bob. Bob is supplied with a 
third measurement in addition to the bases {\l)}i and 
{|e;)};. This measurement is choosen as a third MUB, 

I//) = where = 1 if d is an 

odd prime, and dd = 0 otherwise. 

Alice has (P different encodings but we begin by con¬ 
sidering the d encoding states on the form IV’ooa), where 
we make the ansatz 


IV’ooa) = (|0) + {r + it)\eo) ± (r - it)\fa)) (8) 


where r and t are real parameters and N^^d = 
/ ’\ 1/2 
(l ± ^ ± 2r2 + 2P + 23? (^a) (r^ - P) ± (^„) j 

is the normalization where we have introduced 
= a Z]fc=o This ansatz is made 

for symmetry reasons soon to appear. 

We now claim that given an arbitrary encoding state on 
the form I'i/'ooa), if Bob performs a measurement either in 
the basis {\ei)}i or in the basis {!/;)};, the probability of a 
obtaining the correct value of the corresponding d-levels 
xi and X 2 are always the same. To show this, we only 
need to expand the probability expressions Pi(a:i = 0) 
and P 2 {x 2 = a) given below 


|(eo|l/’OOa)P 


|(/a|l/’00a)P 


^3,d 


^'3,d 


—j= ± r ± it ± (r — it)i^a 

Vd 

^ ± (r ± it)^„ + r-it 


(9) 

( 10 ) 
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where denotes the conjugate of ^a- Observing that the 
real parts inside the modulus squared in both (9) and (10) 
are the same, and that the imaginary parts differ by a 
sign, it follows that the success probabilities of Bob are 
the same, Pi{xi = 0) = ^2(2^2 = a). This symmetry 
property motivates our ansatz in (8). 

We now require that for given encoding state, Bob suc¬ 
cessfully recovers any Xj with the same success probabil¬ 
ity. To enforce this, we compute the success probability 
Po{xo = 0) of Bob. 


a;oa;ia;2 = 00a as 

: So = 0 ^ xo = a (16) 

:xi = 0 ^x 1 =13 (17) 

:S2=a—>'S2=/3 -|-o — (1-1- 5d)a (18) 

where the sums are taken modulo d. It can be proven 
that if Bob measures on the encoding state 

(19) 


Poixo = 0) ^ KOlV^ODa)!^ = ^ (^1 + (11) 

Since we have shown that Pi(si = 0) = P 2 {x 2 = a) for 
all a,d,r,t, we need only to require that ^0(2^0 = 0) = 
Pi{xi = 0). Realizing this constraint using (11) and (9), 
one will find a set of values {r,t{r)) as function of d and 
a such that the state IV'ooa) necessarily fulfills that the 
Bob’s success probability is the same independent of his 
measurement. One is lead to a qudratic equation in t, 
At'^ + Bt + C = 0 with coefficients 


— 1 + 

P-23?(?a) 

(12) 

S = 25 (Ca) 


(13) 


C = r^ [\L\^ + 2^{U + 1 - ^ (3?(Ca) - 1) + ^ 

(14) 


with the familiar solution, now parametrized by r 


t±ir) 


-B ± - iAC 

2A 


(15) 


For given d and a, and arbitrary choice of r (as long 
as the corresponding t± is real), there is a correspond¬ 
ing QRAC with success probability independent of Bob’s 
choice of measurement. However, we are interested in 
QRACs with high success probabilities and it is there¬ 
fore natural to for given d and a, choose the value of r 
that optimizes Bob’s success probability. A key observa¬ 
tion is that the success probability depends on the value 
of the encoding a. In contrast to the case of —)■ I 

QRACs, the average success probability and the worst 
case success probability are not the same for this family 
of -)■ 1 QRACs. 

Before computing the average success probability we 
are required to determine the arbitrary encoding states 
on the form \' 4 ’xoxiX 2 ) which we so far have not consid¬ 
ered. However, general encoding states with good success 
probabilities can be generated from the encoding states 
l^ooa) by acting with the X and Z operators. We choose 
to define the action of these operators on the encoding 


where /(xq, xi,d) = X 1 +X 2 — (l-l-5d)a;o; then the success 
probability is the same as if he measured on \ 1 lJ 00 x 2 ) i-®- 
the success probability is invariant under actions on the 
encoding state lipooa) with X°‘Z^ for a,f3G {0,..., d— 1}. 
Evidently, using actions of X and Z, it is impossible to 
map any |V’oo£C2) to IV’ooa:^) for X 2 ^ ^2 thus we 
are guaranteed to generate d^ distinct encoding states 
with the above method. We can therefore compute the 
average quantum success probability, which by the above 
argument together with the probability expression (11) 
can be written as 



1 

d 


E max 

r,t± 


1 

t±)‘^ 



( 20 ) 


In figure 1 we plot the ratio of the success probability 
of the 3^^') 1 QRACs and the RACs from equation 

1(3). The quantum advantage is maximal for d = 13 at 
which P^d/p^d ~ 1-224 which is signficantly larger than 
the quantum advantage of the known two-level QRAC 
3(2) I fQj. -vvhich p^ 2 Ip 3,2 ~ 1.052. Evidently, the 
quantum advantage can be more than four-doubled by 
considering a high-level QRAC. 

Discussion of optimality. — An important question 
is whether the average success probabilities, p^ d^ 
we have found are optimal, or if there exists other states 
and measurement which lead to better results. To investi¬ 
gate this, we have performed optimization using semidef- 
inite programs (SDPs) [28] , alternating between optimiz¬ 
ing over states and measurements, to find lower bounds 
on p'^ d- We have investigated cases n = 2 (n = 3) for 
d = 2,..., 15 (d = 2,..., 8) for all of which the obtained 
lower bound, up to numerical precision, coincides with 
our presented results. To find upper bounds on we 
have used the intermediate, almost quantum, level of the 
Navascues-Vertesi hierarchy of dimensionally constraint 
quantum correlations [29]. Due to computational limi¬ 
tations, we have only studied the particular cases with 
n = 2 and d = 2,3,4, 5, for which we found that the 
obtained upper bounds coincide with the lower bounds, 
hence proving the optimality of these particular QRACs. 
We also investigated our 3^^^ —1 QRAC by the same 
method. Here, the lower bound, P^d^ 0.6971, is some¬ 
what smaller than the obtained upper bound, 0.6989. 
Evidently, to determine the optimality of this particular 
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QRAC, we need to run higher levels of the Navascues- 
Vertesi hierarchy, which is compuationally demanding. 
The optimality of our 3*^^^ —> 1 QRAC remains an open 
problem. 

Conclusions .— In this letter, we have investigated 
both classical and quantum random access codes with 
d-level communication. We derived the average success 
probability of what we strongly believe is the optimal 
RAC, along with two families of QRACs showing that 
they provide a significant gain in success probability over 
the RACs. Our results show that high-level QRACs en¬ 
able significantly larger advantages over the correspond¬ 
ing RACs than what can be achieved with QRACs using 
two-level systems. The QRACs presented here can be 
applied to a variety of problems in quantum informa¬ 
tion. Nevertheless, this investigation also leaves many 
open questions on the properties of random access codes: 
(1) find an efficient way to generalize our QRACs to ar¬ 
bitrary n. (2) RACs and QRACs can be considered in 
other scenarios e.g. if Bob wants to access a larger subset 
of Alice’s data and/or when Alice is allowed to communi¬ 
cate more than one d-level. (3) Consider scenarios with 
the dimension of Alice’s data set not being equal to the 
dimension of the communicated system i.e. Alice holds 
a string of n d-levels while communicating a single clas¬ 
sical (quantum) d'-level to Bob with d ^ d'. (4) Our 
approach to QRACs is based on MUBs, however it is in¬ 
teresting to ask if this is always an optimal choice of basis 
for any —>■ 1 QRAC, or if there are bases that result 

in QRACs with better performance. 
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